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Abstract 

This talk is an attempt to review all our knowledge on saturation at low x both theoretical 
and experimental, to stimulate a search for saturation effects at THERA. The main goals 
of this presentation are 

1. To discuss an intuitive picture of the deep inelastic scattering that leads to the 
saturation of the parton densities; 

2. To show that the saturation hypothesis has solid theoretical proof; 

3. To report on the theoretical progress that has been made over the past two years in 
high parton density QCD, and on the property of the saturation phase that emerges from 
the theory that has been developed; 

4. To collect all that we know theoretically and experimentally about the saturation 
scale Q s (x) . 



* e-mail: leving@post.tau.ac.il 



Fighting prejudices 



1 Introduction 

At low values of x, QCD evolution, both DGLAPjl] and BFKL|| , predict a striking increase 
of the parton densities which violate unitarity constraints ||. Therefore, interactions between 
partons in the parton cascade, omitted in QCD evolution equations, should become essential to 
slow down the growth of the parton densities. We expect that these interactions will create an 
equilibrium -like system of partons with a definite value for the average transverse momentum, 
which we call a saturation scale (Q s (x)). In other words, we expect a picture of a hadron as 
shown in Fig. 01, |, |, | . 




radius of the target 



Figure 1: Picture of a hadron in the saturation region. 

This talk is an attempt to review all our knowledge on saturation at low x both theoretical 
and experimental, to stimulate a search for saturation effects at THERA. In spite of the fact 
that the high density QCD phase was only briefly discussed in the THERA Contribution to the 
TESLA TDR, everybody knows that if saturation effects are seen at THERA it will be 
the machine, while if not, it will remain one of many. The main goals of this presentation 
are 

• To discuss an intuitive picture of the deep inelastic scattering that leads to the saturation 
of the parton densities; 
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• To show that the saturation hypothesis has solid theoretical proof; 

• To report on the theoretical progress that has been made over the past two years in high 
parton density QCD, and on the property of the saturation phase that emerges from the 
theory that has been developed; 

• To collect all that we know theoretically and experimentally about the saturation scale 
Q s (x). 

Thera are two different ways to reach a high parton density phase: the first, is DIS at low x, 
the second is deep inelastic scattering on a nuclear target in which we have a rather large parton 
density from the beginning, due to a large number of nucleons in a nucleus. The best avenue 
to investigate the high density phase is to use both paths and to measure DIS on nuclei at low 
x. This is one of the THERA options and we will also discuss saturation phenomena for such a 
reaction. 

2 Qualitative Picture of Interaction in DIS at low x 

2.1 Bjorken frame: 

It is well known that deep inelastic scattering is most clearly visualized in a space-time picture in 
the Bjorken frame where the virtual photon has zero energy @. Therefore, in the Bjorken frame 
the electro-magnetic field is a standing wave with wavelength of the order of l/q z for a photon 
with four momentum = (0, q z , 0, 0, ) (see Fig. H). In this frame the fast hadron decays into a 
system of partons. Each parton has a longitudinal momentum 2 = x^P, where Xi is a fraction 
of the energy of incoming hadron carried by the parton, and a transverse momentum Due to 
the uncertainty principle each parton is localized in Azi « l/(xi,P) and, therefore, only partons 
with XiP ~ q z can interact with the photon, since for all other partons the overlap integral is 
very small. In other words, the parton which interacts with the photon has Xi ~ q z /P- Using the 
energy and momentum conservation for the parton - photon interaction one can easily obtainP] 
that XiP = q z /2 which gives Xi = q z /2P = q 2 z /2Pq z = Q 2 /2(P ■ q) = Q 2 /s = XBj at low x. 

The lifetime of the z-th parton is of the order ofr = r r f7 = -^- = where the lifetime in the 

v h ' Pit Pit 

rest frame of the i-parton r r j, = ^ and 7 = The parton that interacts with the photon lives 
for a short time m l/q z and because of this the interaction cannot change the parton distribution, 
and destroys only the coherence of the partons in the incoming hadron. 

Over a long period of time every parton can decay into a large number of partons. Theoretically 
we can only control the emission of the partons with large values of the transverse momenta, as 
the QCD coupling is small for them, and we can safely use the developed methods of perturbative 
QCD. Fig. Q shows the evolution of the partons with definite transverse momenta (p t = Q ^> 1/R) 

1 Energy conservation gives that the energy of the parton i and the recoiled energy are equal (Ei — E'^) , the 
conservation of the longitudinal moment leads to pix = Qz — v'% l- Pi l ~ Vi,h and Pi,L = Qz/2. 
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Figure 2: Parton cascade in the Bjorken frame. 



(with definite size r t = 1/Q -C R) in time ( so called the BFKL evolution). R is the size of the 
hadron. The first stage of the evolution for partons with Xi ~ 1 is not under theoretical control 
and only non-perturbative QCD will be able to give us information on probability (P^(xi,Pit = 
Q -C 1/R ) to find several partons with p ijt = Q and Xi ~ 1 in the hadron. The BFKL evolution 
takes into account the emission of partons with p t = Q ^> 1/R. Fig. g shows that it is natural to 
expect that this emission leads to a considerable increase of the number of partons. The number 
of partons that can interact with the target ( virtual photon) can be written in the form of the 
convolution F BFKL (—,Q 2 ) ® P/(xi,Q). We can obtain the result of the DIS experiment by 
taking the convolution (overlapping integral ) with the photon wave function. In other words, 
the deep inelastic structure function is equal to 

F 2 (x,Q 2 ) =P T {^f,Q) (g) F BFKL (^Q 2 ) P/(x u Q) , (1) 

where P 7 *(— , Q) = |^ r 7 *(— , Q)\ 2 and is the wave function of the virtual photon. 

Recalling that cr(7*, h) = ^ F 2 (xBj, Q 2 ) one can see that the unitarity constraint cr(7*, h) < ttR 2 
leads to a conclusion that the increase of the parton densities due to the BFKL ( or DGLAP) 
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emission should be tamed 0. The simple idea how such taming can occur is clear from Fig. ||]. 
Indeed, if the parton cascade has been measured at early time ( at rather high x ) the density of 
the partons in the transverse plane (see Fig. |2|) is not large, and we have to take into account 
emission, since emission is proportional to the density ( p) of partons ( emission oc p). However, 
the density of partons increase due to emission and at some value of x the system becomes so 
dense that partons cover the hadron disc. In such a situation the interactions of the partons 
start to be essential. These interactions are proportional to the square of the parton density 
since two partons have to meet in one point for this interaction ( annihilation oc («s/Q 2 ) ■ P 2 ? 
where ats/Q 2 is the typical cross section of two parton annihilation in the parton cascade) and 
they cause the number of particles to diminish. Therefore, we expect there to be an equilibrium 
between emission and annihilation in the dense system of partons which can be described by 
simple equation: 

dp _ N c a s ( k bfkl^ p _l^ xp 2) (2) 



d\n(l/x) it \ r Q 2 , 

The first term of this equation gives the BFKL evolution at low x while the second one provides 
a taming of the density increase. Of course, all coefficients in the equation cannot be calculated 
in framework of such oversimplified approach, including numerical coefficient 7. 

The principle prediction of this equation is the saturation of the parton density, namely, the fact 
that the parton density stops increasing. It should be stressed that at any value of Q 2 , even a 
very large value, there exists a small value of x at which we face saturation ( see Fig. ^| ). k oc p 
in Fig. |3| is a packing factor for the partons and we will discuss its value later. 



2.2 Laboratory frame: 

The Bjorken frame is the frame which is best suited for the discussion of DIS in the parton (QCD) 
approach, since both pictures Fig. [| and Fig. |3| give the parton distributions in a hadron, or, 
in other words, term , F BFKL (f:, Q 2 ) (g) P^{x^Q) in Eq. [I]. However, it turns out that some 
properties of the high density parton system are clearer in the laboratory frame where the hadron 
is at rest. As we will see below, the fact that our partons are colour dipoles is easy to demonstrate 
in this frame. The time-space picture of DIS in this frame is shown in Fig. |4j. 

In the lab. frame the fast virtual photon decays into a quark-antiquark pair ( two partons in 
Fig. p. Quark ( antiquark) has transverse momentum larger than Q and exists for a sufficiently 
long time (r « l/ma;)^ If the time r is long enough, quarks (antiquarks) radiate gluons ( as 
shown in Fig. |J) which create a dense parton system in the same way as in the Bjorken frame. 

At first sight pictures Fig. || and Fig. f| look quite different. Of course, the final result of the 
measurement (the total photon-hadron cross section) given by Eq. [TJ , remains the same in both 
frames, but only part of Eq. [I] is shown in Fig. f|, namely, P 7 *(^| i , Q) ® F BFKL (f-, Q 2 ). 

2 The estimates for r we can easily obtain using the uncertainty principle: AE t « 1 where AE is the difference 
in energy between initial and final states. For the virtual photon decay we have AE = <7o — Pi,0+f>2,o ~ qo~qz = mx, 
where pi t o is the energy of produced quark (antiquark). 
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Figure 3: Parton distribution in transverse plane and saturation in the Bjorken frame. 

Therefore, the main difference between these two figures, Fig. and Fig. ^ is the following: Fig. 
D shows all partons which can interact with the photon target, while the system of partons that 
has interacted with the virtual photon is depicted in Fig. [|. One can see that the Bjorken frame 
is much better for describing the parton densities, we will show in the next section the laboratory 
frame is very useful in answering the question: What are these partons in QCD. 

2.3 Colour Dipoles = Partons: 

The advantage of the lab. frame becomes clear if we want to understand how the produced quark 
-antiquark pair( which is a colour dipole) interacts with the target. The main observation is that 
the size (r± in Fig. |^) of the colour dipole or , in other words, the transverse distance between 
the quark and antiquark, is a good degree of freedom, which is preserved by the high energy QCD 
interaction [|, [H| |IT| . Indeed, while the colour dipole is traversing the target, the distance r_|_ 
between the quark and antiquark can vary by an amount Arj_ oc R ^ , where E denotes the 
energy of the dipole in the lab. frame and R is the size of the target (see Fig. || ). Due to the 
uncertainty principle the quark transverse momentum is k± oc Therefore, 
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Figure 4: Parton cascade in the laboratory frame. 



Ar± oc R — ps R r± 

E r± E 



(3) 



r ± 2mE > 2mi?. 

Since r\ ps 1/Q 2 , and recalling the definition of the Bjorken x one can see that 



Ar± 



1 at x 



2mi? 



(4) 



A. Mueller proved two results JO], Q which really showed that the colour dipoles are the correct 
degrees of freedom in QCD at high energies. First, he showed that the gluon structure function 
can be viewed as the interaction of the colour dipole with the target as shown in Fig. |^. Secondly, 
he proved that the BFKL evolution can be rewritten as a decay of one dipole into two dipoles for 
large N c , as one can see in Fig. ||. 

Therefore, we can discuss the high energy (low x) DIS in terms of colour dipoles which interact 
with themselves and with the target. 
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Figure 5: Interaction of a colour dipole with the target (Glauber-Mueller formula) in the laboratory 
frame. 

2.4 Glauber-Mueller formula: 



It turns out that for an interaction with the target we can obtain the simple Glauber- Mueller 
formula which reads f9], |10|, [H], |13| P1 



with opacity 



CTdipole{X,r±) 



n(x,r ± ;b t ) 

1 — e 2 



n 



2 „2 



7T r 



xG 



DGLAP i 



,2 ' 



X) 



(5) 



(6) 



where S{b t ) is the target profile function. In the case of a nucleon target we can use the Gaussian 
form of S(b t ) = (1/iiR 2 ) e~ b ^ R2 while for nuclei we use the Wood-Saxon parameterization for 
S A (b t ). 

3 Giving credit to the authors of Refs. || [h], [l3] we refer to this formula as the Glauber-Mueller formula 
because A. Mueller was the first who proved that the gluon structure function can be described as rescatterings 
of a dipole. This result changed the whole approach to DIS by creating a transparent picture of the interaction in 
QCD at high energies. 
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z 2 2 
x 02 x 12 

Figure 6: BFKL gluon emission as a colour dipole decay. 

Eq. |5| is a solution to the s-channel unitarity constraint for the dipole-target amplitude 

2Ima dipo i e (x,r ± ;b t ) = \a dipole (x,r ± ;b t )\ 2 + G in (x,r ± ;b t ) (7) 
The inelastic cross section is equal to 

<i P oie = J d 2 h (l - e -^W) . (8) 
Opacity Q describes the interaction of one parton shower with the target as one can see in Fig. 

B 

As has been mentioned the real breakthrough was the proof by A. Mueller that the gluon structure 

if 

x G(x, Q 2 



function can be calculated using a similar formula for the colour dipole rescatterings, namely |l I 

8 f 1 dr 2 , 



7T 3 Jx rl 



d%. 1 



(9/4)(l(*,r ± ;b t ) 



(9) 



2.5 Packing factor: 

Eqs. H], § and |^ allow us to introduce a packing factor for colour dipoles in the parton cascade. 
This factor is equal to 

3a s (4-) 7r 2 r 2 / 4 \ 

k = (9/4)0(0^ = 0) = ~^-^t [*G DGLAP ^)) ■ (10) 
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The physical meaning of k is very simple: K = cr dipo i e /7i-R 2 = cr dipo i e (BA)xG(x, Q 2 )/ttR 2 = 
Vdi P oie(.BA) ■ p. Therefore, k is the size of the parton ( or preferable to say its typical cross 
section ) multiplied by the density of gluons in the transverse plane. 

2.6 Observables: 

It should be stressed that all our observables can be calculated if we know the dipole amplitude. 
Indeed, the main obsevables such as the total photon-hadron cross section , the gluon density and 
single diffraction production inclusive cross section, have a very simple relation with the dipole 
amplitude, namely, 

cr(7*p) = J dz J d 2 r ± \^(z,r ± ;Q 2 )\ 2 a d i po i e (x B ,r 2 ± ) ; (11) 

^•« 2 ' = iT^Xf ^,2,1); (12) 
" SD (l"p) = [ d:J ^^(z.rr.Q 2 )] 2 ! d%\a iipa ,,(x,r ± ;b t )\ 2 . (13) 



Eq. |TT|was proven in Refs.|9|, [TO], [IT], [TJ|, while Eq. [12] was first written in Ref.[^JJ and was discussed 

The formula for the total diffractive production in DIS was suggested in 



in detail in Ref.fl3 
Ref. 



The argument 2r\ in Eq. [T2| reflects the fact that actually the rescattering of gluon corresponds 
to rescatterings of two dipoles of the same size, which works effectively as the interaction of one 
dipole but with a size which is a/2 larger (for N c ^> 1). ^ was calculated in Refs.[TT, H|. 



3 Non-linear Evolution 



3.1 The equation: 

Using the colour dipole picture of an interaction and, in particular, the fact that the BFKL 
emission can be viewed as a decay of one colour dipole into two ( see Fig. |] ) we can easily obtain 
the nonlinear equation for the imaginary part of the elastic dipole-target amplitude N(x, r±; bt) = 
Ima dipole (x,r±; b t ). This equation can be written in the form (see also Fig. [7|: 



dN(x 01 ,b t ,y) 2C F 

a S 1 / X 0l\ AT/ 7 a S f i2 X 01 



2 

01 

dy n \ p 2 J 7T J P xg 2 xf 2 

2N(x 02 ,b t ,y) - N(x 2,bt,y)N(x 12 ,b t ,y)) , (14) 



where y = ln(l/x) and we assume that that bt < X02 or/and X\2- 

Eq. p!l has a very simple meaning and actually describes the fact that the dipole of the size xoi 

1 1 2 

decays into two dipoles of sizes X02 and xu with probability |\1/| 2 = r^h as it is shown in Fig. |6]. 



"02 ""12 
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Figure 7: The pictorial form of non-linear evolution equation. 



These two dipoles then interact with the target. They can interact separately and this interaction 
leads to the linear term in Eq. [14]. However, two produced dipoles can interact with the target 
simultaneously generating the non-linear term in the equation. From Fig. |7| one can see that 
this non-linear term takes into account the Glauber correction for the two dipole interaction. The 
minus sign in front of the non-linear term reflects the well known fact that we overestimate the 
value of cross section considering it as a sum of two independent collision, since sometimes one 
dipole happens to be in the shadow of the second one. The linear term in Eq. [14] is the BFKL 
evolution, which describes the evolution of the multiplicity of the fixed size colour dipoles with 
respect to rapidity y. At first sight the linear term sums the leading twist contribution while 
the non-linear one is related to higher twist contributions. However, this is not true. The first 
term ( the BFKL equation ) has also higher twist contributions but with the same anomalous 
dimension as the leading twist ones. On the other hand, as was pointed out by Mueller and Qiu 
the non-linear part contributes mostly to the leading twist. The beauty of the equation is that 
it sums both leading and higher twist contributions in a unique fashion and claims that at any 
fixed Q 2 ( at any short distance ), the higher twist contribution will dominate at sufficiently low 
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X. 



3.2 Brief review of the theoretical approaches: 



Eq. II shows that the problem of high density QCD has been solved from first principles and we 
think that it is instructive to give a brief review of the theoretical approaches that all converge 
to this equation. 

1981 - 1983 GLR pointed out the new phase of QCD-high density QCD, developed picture of 
parton interaction in the Bjorken frame ( see above ), proposed the hypothesis of parton saturation 
and suggested first non-linear equation which sums the "fan" diagrams of Fig. §-a and which is 
actually Eq. |14| in momentum space 0. 

1986 Mueller and Qiu JH proved the GLR equation in the double log approximation of per- 
turbative QCD. 




b) 



Figure 8: "Fan" diagrams (a) and 1/N C corrections to them (b). 



1992 - 1995 J.Bartels[16[ showed that the non-linear equation can be correct only in large 
N c approximation since 1/N C corrections (see Fig. [5|-b ) lead to the interaction of two ladders 
in the "fan" diagrams (see also Ref. |i~7| . Laenen and Levin based on Ref.|1l8| generalized the 
non-linear equation, taking into account l/N c corrections in double log approximation ||19||. It 



turns out that l/A^-approximation works quite well in this problem and can be treated using the 
generating function formalism. 

1994 L. McLerran and R. Venugopalan || noticed that at high density the gluonic fields are 
strong ( G^y ~ 1/g where a s = g 2 j^) and, therefore, one can approach the high density QCD 
using the semiclassical gluon fields. Based on the space-time structure of the parton cascade at 
low x, they built the effective Lagrangian for high density QCD. 
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Figure 9: The space-time structure of the parton cascade at low x. 



Indeed, the main and very important feature of the parton cascade, shown in Fig. [| is the fact 
that a parton with higher energy in the cascade lives much longer than a parton with lower 
energy. We follow the parton emitted at time t' 2 . This parton lives a much shorter time than 
all partons emitted before. Therefore, all these partons will have enough time to form a current 
which depends on their density. Since the density is large we expect that the current is a classical 
current. Finally, the Lagrangian of the interaction for the parton emitted at time t' 2 can be written 
as 

L(p) + j„ ■ + L(A) , 

where A is the field of a parton emitted at t' 2 . However, we can consider a parton emitted at 
t = t' 3 and include the previous one in the system with density p. The form of the Lagrangian 
should be the same. This is a strong condition (equation) on the form of the effective Lagrangian, 
so called Wilson renormalization group approach. This is the beautiful idea of the McLerran and 
Venugopalan which leads to Eq. [14], as we will show below. 

1996 I. Balitsky |20| proved the non-linear equation in Wilson Loop Operator Expansion at 
high energies developed by him. Unfortunately, his paper was not noticed by the experts in 
the field including me. It should be stressed that he also gave an operator proof of the BFKL 
equation. 

1997 Ayla, Ducati and Levin suggested non-linear equation [I3| which differs from that of 
Eq. [14]. They used the double log approximation but summed two DLA contributions 

((a s ln(l/x) ln(Q 2 /A 2 )) n and (a s ln(l/x) \n(Ql(x) / Q 2 )) n . It turns out that their equation is just 
the same equation as Eq. [13] but written for the opacity fl instead of N = (1 — exp(— ^)). AGL 
determined all numerical coefficients and discussed the initial condition of Eq. [14] which we will 
consider later. 
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1999 -2000 Yu. Kovchegov[21] proved Eq. 14 in the colour dipole approach [12]. It should 



be stressed that he not only gave the derivation which we have discussed, but he found a correct 
observable which enters the equation (N) and he suggested an initial condition that we will discuss 
below. 



2000 M. Braun |22] calculated the "fan" diagrams of Fig. [8|-a in the BFK1 kinematic region, 
using the triple ladder vertex of Refs. P3fl . 

2001 Iancu, A. Leonidov and L. McLerran [Q (see also Ref. |[25|| ) proved Eq. [H]m the effective 
Lagrangian approach. Their proof is based on long and successive development of the effective 
Lagrangian approach exploited in Refs. [26] 



3.3 Initial conditions: 

One can see that Eq. |TJ] does not depend on the target and the dependence on the target comes 
from the initial conditions at some initial value of x = Xq. For a target nucleus it was argued 
T3, ^T[ that the initial conditions should be taken in the Glauber - Mueller form (see Eq. |^), 



namely, 

/-~* 71 f ^(^01 ' X ~ X '^t ) 

N(x 01 ,x = x ,b t ) = N (x i,x = x ,b t ) = 1 - e a . (15) 
The value of Xq is chosen in the interval 



i 

Q s < Xq < 



2mR 



where R is the radius of the target. In this region the value of Xq is small enough to use the low x 
approximation, but the production of the gluons (color dipoles) is still suppressed as as ln(l/x) < 
1. Therefore, in this region we have the instantaneous exchange of the classical gluon fields. Hence, 
an incoming color dipole interacts separately with each nucleon in a nucleus (see Ref. ||27|| ). 

For the hadron we have no proof that Glauber- Mueller formula is correct. As far as we understand 
the only criteria in this problem (at the moment) is the correct description of the experimental 
data. We described PS| almost all available HERA data using Eq. ||, and we feel confident using 
Eq. |5| as the initial condition for Eq. |14[ It should be stressed that the experimental data on 
dF 2 /d In Q 2 provides direct information on the integral over b t for N since 

dF 2 (x,Q 2 ) Q 2 



d\nQ 2 3tt 3 



J d 2 b t N(xA/Q 2 M) + 0(l/lnQ 2 ) . (16) 



Choosing x — we can make the initial condition practically independent of the parameter- 
ization of the xG DGLAP since all available parameterizations give the same prediction, even for 
the gluon density in this x - range. 



3.4 Theory status of the approach: 
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3.4.1 Parameters of the approach: 



The master equation (Eq. 0) as well as the Glauber - Mueller formula (see Eq. | and Fig. || ) 
sums all diagrams of the order of 

U|(-) A j with A oc a s (17) 

It means that starting from as ln(l/a;) ~ ln(l/as) the corrections related to the rescatterings 
of the partons and their interactions become essential. On the other hand, the next to leading 
order corrections to DGLAP or BFKL equations produce a value of A in Eq. [Tj] calculated 
up to a s accuracy ( A = C\a$ + C 2 a s ). Such corrections start to become important only for 
as ln(l/x) > l/«5 or , in other words, the next to leading order corrections have to be calculated 
after the rescatterings and parton interactions, included in the master equation, have been taken 
into account. 

Therefore, the calculations of the next to leading order corrections for the linear 
evolution equations can be considered as a QCD motivated model in which crucial 
ad hoc assumptions have been made, namely, that the parton interactions and their rescatterings 
are small numerically. 



3.4.2 Large N c approach: 

As has been discussed the master equation (Eq. |14]) has only been proved in the large N c limit 
of QCD. The key question arises: could we develop a selfconsistent large N c approach for the 
DIS at low x. For JV C > 1 we assume that &s = ds^c ~ 1 while as <C 1, we assume all 
contributions that are proportional to as to be small so that we can neglect them. At first sight 
we have a problem with developing such an approach for high energy scattering. Indeed, the first 
Born diagram ( see Fig. |TD| ) turns out to be of the order of (Nf/N c )a s [] where Nf is the number 
of the quarks inside the proton. We can consider this parameter as being small in the large N c 
limit. However, it is well known that a nucleon consists of the N c ( Nf = N c ) number of quarks. 
This observation makes (Nf /N c )a s a 'small parameter' of the order of a s - 

The amplitude A(x 01 ,6 t = 0,y) « a s /'KR 2 where R is the radius of the proton for N c 3> 1. 
For large N c the interaction of two quarks is proportional to as and can be considered as small, 
but in the Hartree-Fock approximation it leads to a potential proportional to asNf = as N c . 
Therefore, R cx as- 

Finally, we conclude that the master equation sums all contributions of the order (as ln(l/x)) n 
and gives a selfconsistent large N c limit for the DIS with the nucleons. In such an approach the 
DIS with mesons is suppressed by extra factor 1/N C . However, the equation itself depends only 
on as, and 1/N C suppression for the mesons means that the initial condition could be written 

4 This estimate comes from a| (j) A ~ 1. 

5 In large N c limit we can look at gluon as the quark-antiquark pair and Fig. [1^ shows only one quark loop, 
which gives factor 7V C , for the diagram which is of the order of a|. 
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Figure 10: The Borm approximation for the interaction of the virtual photon with the proton in large 
N c limit of perturbative QCD. 



only as the first term of the Glauber- Mueller formula (see Eq. ^ ). Using this initial condition, 
we find the first 1/N C corrections to the amplitude which will give the main contribution for the 
high energy ( low x ) asymptotic. 

It should be stressed that we still do not need to take 1/N C corrections related to more complicated 
diagrams of the Fig. [5|-b type into account. We would like to emphasize that the regular procedure, 
of how to calculate 1/N C corrections, has been developed |L9], [25[ and the first estimates have 



been made [19 . 



3.4.3 The region of applicability: 

To evaluate the value of energy up to which we can use the master equation as a good theoretical 
tool for finding the high energy asymptotic, we need to estimate the contributions that have been 
neglected. We discuss all of them of the order of 1/N C , and they can be divided in two major 
classes: 

• The first one is the interaction between two "ladders" of type Fig. |]-b . They are pro- 
portional to as/N^ ln(l/x) ||16|| , and therefore, we have a first restriction on the value of 

X' 

-| ln(l/x) < 1 or ln(l/x) < ; (18) 



• The second constraint comes from so called enhanced diagrams (see for an example Fig. [IT 
). The ratio of the contribution of this diagram to the "fan" diagram of Fig. |8|-a is pro- 
portional to as /N% ln(l/x) which leads to the same value at the highest energy, as does 
Eq. [LJ. 

Eq. [18] gives a sufficienly large value of energy upto which we can trust the master equation. The 
actual situation is much better. Indeed, the nonlinear equation leads to N{^x 01 ,b t ,y) — > 1 at 
high energy, and we will argue below that this asymptotic limit will be reached at much lower 
energy than one denned by Eq. [TR Since iV(xoi, b t , y) = 1 is the unitarity boundary for the 
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Figure 11: The first enhanced diagram in the parton cascade. 



scattering amplitude, the 1/N C corrections cannot modify this result but could lead to a slightly 
different value of the saturation scale. 

It should be stressed that the regular procedure of how to include l/N c corrections to the master 
equation, has been developed [20|, |24|, |25|] and even first estimates of the quantative effect of 
the 1/N C corrections have been performed fTTfl . 



4 Saturation scale Q s (x) 



4.1 Simple estimates: 

The simplest estimates for the saturation scale come from the equation that the packing factor k 
is of the order of 1 [11, 13, p5| . Fig. 12 shows the solution of the equation 

k(x,Q 2 s (x)) = 1 (19) 

As one can see Q s can be rather large even in HERA region. However, this estimate depends 
crucially on the value of xG DGLAP . 



4.2 Analytic solution of the non-linear equation: 

Eq. |14] has been solved analytically with the simplified version of the BFKL kernel |3(], [H], p5| , 
namely, the Mellin transform x(f) was taken as j + jzj instead of x(f) = 2^(1) — ^(/) — ^(1 — /)> 
where ip is the derivative of the log of Euler gamma function. In terms of s-channel resummation 
Eq. [14] with this kernel sums two kinds of double logs: (asln(l/x) ln(Q 2 /A 2 ))" for Q 2 > Q 2 s (x) 
and (a 5 ln(l/a;) hi(Q 2 s (x) / Q 2 )) n for Q 2 < Q 2 (x). For large Q 2 Eq. |T3|can be reduced to the linear 
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equation and can be solved using the trajectory (characteristics ) method [||, [HJ. The structure 
of the trajectories both for linear and nonlinear equation is shown in Fig. 13. 

The idea is to interpret the last linear trajectory ( line with u cr = 2 as in Fig. |l3|-b ) of the linear 
equation which can be treated without nonlinear corrections, as the saturation scale. Doing so, 
we obtain RO] the following expression for the saturation scale 



Qi 



Ql 



[ X 



x ■ \ — 



Xq 
X 



4N c as 



(20) 



where Q 2 s {x 



Xq)) is a saturation scale in our initial condition. 
Therefore, we expect power-like rise of the saturation scale at low x. 



4.3 Phenomenological saturation scale: 

Golec-Biernat and Wuesthoff [|33] suggested one could extract the value of Q s (x) from HERA data 
assuming that the saturation region has been reached at HERA. Surprisingly they managed to 
describe almost all HERA data using a simple parameterization for the saturation scale, namely, 

Ql{x) = (IGeV*) • (^) A (21) 

with A = 0.288 and x = 3.04 x 10~ 4 . 
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Trajectories for linear DGLAP eq. y=( %? r + Vo q2 = 9§ y ) 

a) b) 
Figure 13: Trajectories for linear (a) and non-linear (b) equations. 

4.4 Saturation scale from numerical solution of the non-linear equa- 
tion: 



In Ref. []34| an attempt was made to solve the non-linear equation numerically (see also Ref. |22| 
), starting from initial x = xo = 1CT 2 . Denning the saturation scale as a value of r\ = A/Q 2 s (x) 
at which the imaginary part of the elastic amplitude for the dipole-target scattering is equal to 
1/2 

N(r ± = 2/Q s ,x) = ^ 
the saturation scale shown in Fig. [14| was calculated. 

One can see that the saturation scale reaches a sufficiently large value of the order of Q s (x) ~ 
14GeV at x = 10~ 7 . It should be stressed that even at x — 10~ 5 which is in the typical range 
for THERA, the saturation scale is approximately 2 v 3 times larger than the Golec-Biernat and 
Wuesthoff estimates. 



4.5 Saturation scale from DIS with nuclei: 



The first estimates for the saturation scale in DIS on nuclei, presented in Fig. [T5|, shows that the 
nucleus target gives a promising avenue to increase the parton density without requiring a region 

-3 



of extremely low x. Comparing Fig. [L5| with Fig. [14] we see that Q s (x) at x 
almost twice larger than for a nucleon. 



10 for gold is 
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Figure 14: The satu- 
ration scale 
Q s (x) plotted 
as a function 
of Igx = 
log w {x). GW 
denotes the 
saturation 
scale from 
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A new scaling in the saturation region. 



The simple picture of the hadron in the saturation region shown in Fig. [I], leads to a new scaling 
phenomenon in this saturation region. We expect that the parton densities as well as cross 
sections are not functions of two variables:x and Q 2 , but they depend only on ratio Q 2 /Q 2 (x) 
H [35], H [15], 3~0fl . Stasto, Golec-Biernat and Kwiecinski J36| found that this scaling is valid for 



HERA data at x < 0.01. Fig. [16] summarizes the situation and, here, I would like to comment on 
the theoretical arguments for such a new scaling scheme. 



5.1 Simple arguments for a new scaling: 

Let us rewrite Eq. [14] in momentum space where it looks simpler 
dN(Q,y,b t 



dy 



- = «S (/ K(Q, Q')N(Q>, y; b t ) - N 2 (Q, y; b t 



For moments N(u, Q) = N(u, Qo) e^ )lnQ ' 2 Eq. [22] can be reduced to the form 

ue l{u>)lnQi = ^^HZnQ^ _ N ^ J ^ e ( 7 (a/) +7 (u,-a/)foQ 2 

The integral over u' in Eq. ^3] can be evaluated by saddle point method 



(22) 



(23) 



Therefore, we can see two different regions in the solution to Eq. large Q 2 , where the non- 
linear term is small and can be neglected, and low Q 2 where linear and nonlinear terms should 
be of the same order. It gives |35| 



7 (w) = 2 7 (w/2) 



j(uj) = Coj 



(24) 
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Figure 15: Saturation scale for DIS with nuclei: A (a) and x (b) dependencies 



Using Eq. 24 we obtain for N(y, Q 2 ) 

N(y,Q 2 ) = J due-^ + ^ )laQ2 = J duoe^- y+ChlQ2) , 

which means that N(y, Q 2 ) depends only on one variable |35| 

-y + C\nQ 2 = \n(Q 2 s (x)/Q 2 )=Z. 
It turns out that C = 41 ^ as in double log approximation P^fl . 

5.2 Scaling solution: 



(25) 



(26) 



Eq. 26 gives the equation for the critical line Q 2 = Q 2 {x) in Fig. [3]. A more refined approach to the 
solution of the master non-linear equation, allows us to understand how far we have to move from 
the critical line to see this new scaling. It turns out that the scaling is not applicable |3(| only in a 



narrow (at low x ) band along the critical line with the width \n(Q 2 (x)/Q ) ~ l/4ag ln(l/x) 1 
for as ln(l/x) 3> 1. 

It should be stressed that two quite different approaches: one is the analytic solution of the 
master equation |$0| and the second is the solution of Wilson renormalization group equation for 
generating function ||25|| , lead to the same answer and the same picture for the new scaling. The 
scaling solution is given in Fig. the difference between N(z, b t = 0) and a is due to the 

integration over impact parameter bt. 
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5.3 Scaling violation: 

The question which we wish to address is how large is the scaling violation for Q 2 > Q 2 s (x). In 
Ref.[pO[] one can find first estimates of this scaling violation (see Fig. [18] 

One can see that at all reasonable values of y = la(l/x) and £ = \n(Q 2 (x)/Q 2 ) the violation is 
less than 30%. It means that if we want to have a tool to measure Q 2 {x) as a value of Q 2 at which 
we see a scaling violation, we have to establish a scale to an accuracy of less than 10%. We think 



that it is instructive to examine Fig. [19| from this point of view. One can see from this figure that 
the accuracy of new scaling phenomena in HERA data is low, and we cannot use these data to 
measure the value of the saturation scale. 



6 Parton densities at THERA and LHC energies 
6.1 Analytic estimates: 

All three theoretical analyses || |30], |25j of the master equation give the same result for the dipole- 
target amplitude N(x,r±;bt), namely in the region of low x N(x,r±;bt) — > 1 in accordance 
with the unitarity constraint. Eq. [1^ leads to gluon structure function which can be calculated 
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Figure 17: (a) Dipole - target scattering amplitude N(z) at bt = and (b) dipole - target cross 
section cr(z') in the scaling approximation versus scaling variable (a) z = ln(r^ Q 2 (x)) and 
(b) z' = ln(Q 2 (x)/Q 2 ): Solid line is the scaling solution of the master equation, dashed 
line is a Golec-Biernat - Wusthoff model as explained in text and dotted line is the z>l 
asymptotic calculated in the first paper of RefpOfl . 



as follows 



xG(x,Q 2 ) 



Ti" 



1 dx' [°° drj_ 



x 



4/Q 2 rl 



d 2 b t N(x' } r ± ;b t ) 



— ~ 2 Q 2 R 2 rdy', (27) 

7T 2 Jy 

where y cr is the solution of the equation Q 2 (y cr ) = Q 2 and R is the size of the target. Therefore, 



the answer for the asymptotic is clear, and it can be rewritten in the form 0, RTj 



xG(x, Q z 



N 2 



1 1 



47r N c as 



— R Q HQ 2 s (x)/Q 2 



(2* 



In Eq. ^8| we assumed that the b t -distribution does not depend on x. Actually this is not correct 
and we have some shrinkage of the diffraction peak |[30|| . All numerical coefficients in Eq. ^ are 
fixed in the double log approximation and have to be checked to a better accuracy. 



6.2 Numerical results: 

The analytic approach leads to an understanding, but also gives a certain check of a numerical 
procedure for solving the master equation. At the moment, we have two attempts to solve the 
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Figure 18: The ratio 5N(y,£)/N(£). 

equation numerically |22|, Q which can illustrate the effect that the non-linear evolution will have 
for the extrapolation of HERA data to higher energies ( lower x). (see Fig. 

We can learn at least two lessons: first, the approximate models cannot be used for the predictions 
at higher energies; and second, the saturation phenomenon could be rather strong both at THERA 
and LHC energies. Such a strong modification of the r± - behaviour due to non-linear corrections 
reflects in Q 2 and x behaviour of gluon density as one can see in Fig. |21 

Fig. £H] shows that the taming of the parton densities growth included in the master equation 
can easily reduce the value of the gluon density at THERA and LHC energies by 24- 3 times, at 
different values of Q 2 . 



7 Summary 

7.1 From the first principles: 

The main message that I want to deliver in this presentation is very simple: Due to the hard work 
of number of experts, the theory of high parton density QCD has been established, and this theory 
includes the proof of : 
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Geometric scaling 




Figure 19: The HERA data, 
plotted function 
of t = Q 2 /Q 2 s (x). The 
solid curve is the pre- 
diction of the Wusthoff 
and Golec-Biernat 
saturation model, 
which was calculated 
by K. Golec-Biernat 
specially for THERA 
WS. The two dashed 
curves give the pre- 
diction of the DGLAP 
evolution (parameter- 
ization GRV'98) for 
the THERA kinematic 
region. 



• The nonlinear evolution equation for the dipole-target amplitude at fixed b t (Im a(x, r±;b t ) 
N{x,rr,b t )); 

• The saturation of the parton densities as x — > which means that N(x, r_|_; b t ) — > 1 at 
low x; 

• The new scaling phenomena N(x,r±;b t ) = N(rj_ ■ Q 2 (x),b t ) for Q 2 < Q 2 (x); 

• The sharp increase of the saturation scale Q 2 (x) in the region of low x; 

• The importance of the saturation effects in taming the growth of the gluon density at 
THERA and LHC energies. 

We would like to stress that during past two decades we have developed the perturbative QCD 
methods [21], based on the correct degrees of freedom: colour dipoles at high energies []TT|1 and 



created new operator methods of tackling this problem, such as the Wilson Loop Operator 
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Figure 20: The comparison between the solutions N, Glauber- Mueller formula [28], and Golec-Biernat 
and Wuesthoff (GW) model. The four curves correspond to two different solutions N as= Q^5 
(large dashes), N as - r unning (small dashes), Nqm (continuous line), and Nqw of the GW 
model (dots). 



Expansion ||20|| and the effective Lagrangian approach^, [24], [2~5f . My personal feeling is that the 
future is in the operator methods since they provide a possibility to describe the matching of the 
high density QCD with real non-perturbative QCD with large coupling constant in a unique way. 
However, I would like to emphasize the positive aspects of the pQCD approach : a clear physical 
picture in the pQCD calculations and, because of this picture, the transparent understanding of 
the meanings of observables which we calculate in pQCD. 

The whole development of this new area of theory I consider as a triumph of approach of my 
teacher, Prof. Gribov, to high energy physics, which could be formulated as "Physical picture 
- first, mathematics - after if ever"|38|]. Indeed, I am very certain that the remarkable 
progress, that we see now, was possible only, because the picture that has been discussed in the 
first section, was correct from the beginning. 



7.2 Our prejudices: 

I decided to finish my paper listing the prejudices (some of them ) that I have fought during the 
two decades: 

• The DGLAP evolution equation is more fundamental and has better proof than non-linear 
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Figure 21: The function xG is plotted versus lg(x). The small dashes correspond to the solution with 
the running as, while the large dashes are used for the constant as = 0.25. The continuous 
line is the GRV parameterization. 

one. It is not true because 

— The GLAP equation yields the parton densities which violate the unitarity con- 
straints; 

— The next to leading order corrections to the DGLAP evolution equation have to be 
included at energies ln(l/a;) > ( see section 3.4.1) which are much higher that 
the energy ln(l/x) > ln(l/ as) /as when the parton rescatterings and recombinations 
start to be essential; 

— It does not contain any proof of why the higher twist contribution is smaller, than 
the leading twist one, which has been incorporated in these equations; 

— These equations demand initial conditions which we do not know how to treat the- 
oretically or extract experimentally. 

• By choosing the initial value for the DGLAP evolution Qq to be large enough, we can claim 
that the higher twists contributions are small. It was shown that the anomalous dimension 
of the higher twists is much larger that the leading one [|16| , |i~7l and at any value of Q 2 at 
low x, the higher twist will exceed the leading twist (DGLAP) contribution. Therefore, the 
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practical procedure of solving the DGLAP evolution equation, cannot be considered to be 
free of justified criticism ; 

• The calculation of higher order corrections of pQCD will improve our description of 
the experimental data. The pQCD series are the asymptotic series and they have intrinsic 
accuracy, which cannot be improved by calculating additional orders in as, and this accuracy 
is rather poor in the region of low x. To evaluate errors we have to use the next order 
calculation in pQCD as an error, but with the same initial conditions. 

7.3 The map of QCD 



npQCD Saturation Scale 




Colour 
Glass 

Condensate 



A 



hdQCD 



pQCD 
BFKL 



DGLAP 



Colour 

Dipole 
Gas 



CD 



1 



Figure 22: The map of QCD. a s 
is the running QCD 
coupling constant, 
k = (3ir 2 a s /2Q 2 ) ■ 
{xG{x,Q 2 )/ttR 2 ) 
is the packing fac- 
tor for the gluons. 
The saturation scale 

Q 2 = Qs( x ) 1S tne nne 

on which k = 1. 



'/oc s 
1/k 



Fig. |22| displays the QCD map which is as the result of this approach. Three separate regions 
are distinguished by the size of the variables Q 2 and x. Each region corresponds to quit different 
physics. 

1. Perturbative QCD domain where the constituents are of small size and their density is 
rather small so that we can neglect their interactions (recombinations) . The packing factor k < 1 
, and we can view this region as the traditional region for the linear QCD evolution equations ( 
the DGLAP and BFKL ones ). In this region the parton system appears as a Colour Dipole Gas, 
since we have here a rather dilute system of colour dipoles. 

2. High parton density QCD domain in which the constituents are still small and we can use 
weak coupling methods, but their density is so large that the gluon field in this region becomes 



strong. We call the parton system in this domain a colour glass condensate |37|], because gluons 
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have colour, their fields evolve very slowly relative to the natural scale and are disordered. In 
simple words glass is a solid on a short time scale and a liquid on long time scale. In lectures by 
L. MacLerran |37|]the reader can find the arguments why the parton system in this domain looks 



like a glass. We have a condensate because of high parton density and strong gluon fields in this 
region. 

3. Non perturbative QCD domain in which the QCD coupling is large (a s ~ 1). In this 
region the confinement of quarks and gluons occurs and a real non perturbative approach should 
be developed here. The key words for this domains are QCD vacuum, lattice QCD, effective 
Lagrangian approach, QCD - string theory duality and other lofty words about possible solutions 
of non perturbative QCD. 
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